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Abstract 

We consider the Helmholtz equation in the half space and suggest two meth- 
ods for determining the boundary impedance from knowledge of the far field 
pattern of the time-harmonic incident wave. We introduce a potential for which 
the far field patterns in specially selected directions represent its Fourier coef- 
ficients. The boundary impedance is then calculated from the potential by an 
explicit formula or from the WKB approximation. Numerical examples are given 
to demonstrate efficiency of the approaches. We also discuss the validity of the 
WKB approximation in determining the impedance of an obstacle. 

1 Introduction 

Various impurities such as gases, dust, cracks, etc., on the surface of a body subject 
to an incident wave can be modeled by the impedance boundary condition [1]. The 
detection of these inhomogeneities using nondestructive testing is then reduced to the 
reconstruction of the impedance from the measurements of scattering field [2] . Optical 
scanning of the surface of silicon wafers used for quality control in semiconductor 
industry [3] is one of the possible applications of this method. 

We consider the scattering of an incident time-harmonic plane wave from the bound- 
ary of the half-space M.^ = {x = (xi, X2, x^), x-^ > 0}. The problem is described by the 
Helmholtz equation 

- An = k'^u, X3 > 0, (1) 
with the impedance boundary condition 

+ ik-f{x')u\^^^Q = 0, x' = {xi,X2,0), (2) 
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where 7(0;') is the surface impedance with a bounded support supp7 C [—1, 1] x [—1, 1] 
and u = u{x) is the superposition of the incident, reflected, and scattered waves 

u{x) = e''''^ + + ^{x). (3) 

Here k = (/ci, /c2, /cs) is a vector such that \k\ = k, k* = {ki, k2, —k^), and function 
^l>{x) satisfies the radiation condition 



i/j^x) 



iklxl 



\X\ 



f{k,x)+0 



x = -^eS\ (4) 



\x\ 



The inverse scattering problem for ([I])- ([2D consists in determining the impedance 7(0;') 
by the far field pattern f = f{k, x) when k is fixed and S G S^. 

In the next section we introduce a modified potential v and express the impedance 
7 through V using an explicit formula. The mapping v ^ f is linear. Hence, the 
initial nonlinear inverse problem is split into two steps: solution of a linear problem 
(restoring v from /) and application the explicit formula. The similar approach was 
used in the discrete counterpart of the problem [1]. This approach does not formally 
require k ^ 1. We also modify it for large k using the WKB method. 

In the case of a bounded obstacle, the WKB method allows one to connect the 
impedance 7 with the asymptotic expansion of the far field pattern (see [S]). We per- 
form a simple numerical calculations in order to find the range of parameters for which 
the WKB method can be used to determine the impedance. The inverse impedance 
problem has been considered in [6]- [7] for general bounded obstacles. Our assumptions 
simplify the problem, and as a result its analytical and numerical solutions become 
easier. 

Note the difference in WKB approach in the inverse impedance problem for the 
half space and a bounded obstacle. In the latter case the asymptotic expansion of the 
far field in a given direction is determined by the value of the impedance in a specific 
point if the obstacle is convex. This is not true for the half space. 



2 Explicit formula for the impedance 

We reduce the problem to the whole space M.^ by extending function u evenly through 
the boundary X3 = for 0:3 < 0. Then equations ([I])-([2]) are replaced by the Schrodinger 
equation 

{-A + q)u = k'^u, xeR^, (5) 

where potential q{x) = —2ik'-f{x')6{xs), and 6{x) is the Dirac delta-function. 

Substituting ([3]) into ([5]), we obtain that the scattering solution ip{x) satisfies the 
equation 

(-A + q{x) - = -q{x) (e*'^-" + e''"*'") = -2q{x)e"''-'^' . (6) 

Equation ([6]) is uniquely solvable if ip satisfies the radiation conditions (jlj). From ([6]) 
it follows 

(-A - P)^ = -q(x) (2e"'' '^' + V) = 2ik-f{x') (2e"'' '^' + ^P{x')'^ S{x3). (7) 
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Let us denote by c{x') the coefficient of S{x') in the right hand side of ([7]) 

c{x') = 2ik^{x') (2e'^'-''' + ^{x')^ . (8) 

In this notation equation ([7]) has the form 

(-A - k^)^ = c{x')5ix-i). (9) 

Observe that coefficient c{x') vanishes outside the support of '~f{x') and hence sohition 
of equation ([9]) can be written as 

Hx)= I G{x-y')c{y')dy\ (10) 

^ik\x-y\ 

where Gix — y) = —-, r is the Green's function of Form (fTOil and (El) we 

4:71 \x — y\ 

obtain equation for determining c{x') 

c{x') + q{x') I G{x' - y')c{y') dy' = -2g(x')e^^'-^'. (11) 

J supp 7 

Finally, it is convenient to introduce a modified potential v{x') as 

v{x') = -c{x')e~''''''' . (12) 

TT 

Then (fTTl) becomes 



v{x')2ik-^{x') / G{x' - yy'' ■'^y'-'^My') dy' = — lix'). (13) 

J supp 7 ^ 

Thus, if f (a;') is known, one can find 7(2;) from (fT3l) from the formula 

/\ iv(x')k~^ . 

-^(^ ) = -4.- + 2 G(=.' - ^')e-'-(-'--).(2/') rfy" "^^"PP^- ^^4) 



In the next section, we will describe a method of determining v{x') from the far 
field pattern f{k,x) with fixed k, and this will complete the solution of the inverse 
impedance problem. 

3 Calculation of the modified potential 

Equation flTU]) contains Green's function of a shifted argument whose asymptotic be- 
havior has the form 



k\x\ 



1 p«« 

G{x-y) = -^ 
An \x\ 





1 - 


-(r)1 













\x\ — > 00. (15) 
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Substituting it into (jl]) and ffTOj) . we obtain the following representation for the far field 
pattern f{k, x) 

f{k,x) = \ I e-'^'^'-^yy' v{y') dy\ (16) 

J supp 7 

where y' = (|/i,?/2,0) and x = x/\x\. Our next goal is to select directions x so that 
the integral 0161) would represent the Fourier coefficients of function v{y'). 

To this end, we write down the incident vector as = /c(cos y?!, cos ip2) cos ^93), while 
X = (cos6'i, cos6'2, cos6'3). Then ( TTB]) becomes 

f{k,x) = ^J j e-^*^[^('=°'^i-™'^i)+J'(™'^2-™'^2)li;(x,|/)rfx(i|/. (17) 

Expression (fT7|) can be associated with the Fourier coefficients of v{x,y) if angles 
6'i = 6*1,711 and 62.n2 are chosen in such a way that 

A;(cos6'i„i — cosy^i) = nrii, (18) 
/c(cos6'2,n2 - COSV92) = vrn2, (19) 

where rii, n2 = 0, ±1, ±2, . . ., and 

k k 

{1 + cos (fi) ^ rii ^ — {1 — cos (fi) , i = 1,2. (20) 

For those directions defined by the angles 6*1 and 6'2,n2 ? the measured far field pattern 
fni,n2 will be the Fourier coefficient in the expansion of the modified potential v{x,y) 

fnun, = \ f ^ e--^("^^+"^^) v{x, y) dxdy. (21) 
Hence, v{x,y) has the following Fourier series representation 

v{x,y)= ^ /„,.2 6'^^^'^^^+"^^^ (22) 

ni,n2 

Formula (I14p along with (122]) provides the solution of the inverse impedance problem. 



4 Asymptotic solution 

Now we are going to modify the previous approach assuming A; ^ 1 and using the 
WKB approximation. The scattered wave ip{x) satisfies the Helmholtz equation in the 
half space M.^ = {x = (xi,X2,a;3), X3 > 0} 

- Ai> = k^ij, X3 > 0, (23) 

and the boundary condition 

^^^3 + ikj{x')ip{x') = -2ikj{x') e''''-'^', x' = (xi, X2, 0). (24) 
13=0 
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In order to find asymptotic behavior of tp{x) for large k, we will use the WKB ap- 
proximation of ip{x) in a neighborhood of support of 7(2;'). We will be looking for 
expansion of ip in the form 

ij{x) = e**=* =^ Yl *"(^) (^^)~" ' ^* = (^1' ^^2, -^3). (25) 

n=0 

Coefficients in this expansion can be found explicitly. Substituting (!25|) into equation 
( !23|) and equating the coefficients of like powers of k, we obtain a recurrence system of 
differential equation for "^ni^) 

p.V^o = 0; (26) 

2fc^ ■ + A^„_i = 0, n^l, (27) 

where k* denotes the unit vector in the direction of vector k*. From the boundary 
condition ( 124|) . one can find the initial condition for \E'n(a3) and thus determine all the 
coefficients \E'„(a;). In particular, ( 126|) implies 

*o(a;i, xs) = $0 (xi + ^2 + ^) , (28) 



where $0 is an arbitrary differentiable function. From flMj) it follows that 

,T, f a\ 27(xi,X2) 

'^°("^'"^'°^ = - 7(xi,x.)-fc3fc-^ ' (29) 

and hence 

/ X 2-f{xi + kik^^X3,X2 + k2k^^X3) , . 

^0[Xi,X2,X3) = ; ^ — j ^ — j -. (30) 

7(xi + kik^ X3, X2 + k2k^ X3) - k^k ^ 

Thus, using expansion fl25l) and relations ([8]) and f|T2l) . we obtain the following asymp- 
totic representation for the scattering amplitude / fITBl) through the boundary impedance 

7 

fik, x) = -^ I e-(^-^)-^' . ^ ^ dy' + O {k-') . (31) 

Jsupp 7 7(y ) -hk ^ 

If we choose the direction of measurements x of the far field pattern the same as before 
in f|T8|) - fl20|) . then the value f{k,x) becomes proportional to the Fourier coefficient of 

^(^-ksk-')-' 

fn.,n, = e-'^^(-^^+-^^) , 1 dy^dy^ + O {k~') ■ (32) 

J-iJ-i 7(^1,2/2) - k^k-^ 

Applying the inverse Fourier transform, we obtain 

-f{xi,X2) ni 



7(xi,X2) - ksk 1 Aks 



T = jrYl e-^(™^'^+"^^) + O {k-') . (33) 



This equation can be solved for 7. Hence the boundary impedance can be restored 
using the values of the far field pattern in the specific directions given by (fT8l) -( |T9l) . 
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5 Scattering from sphere 

In the case of convex body Q, there is a direct asymptotic relation between the far field 
pattern and the boundary impedance for large values of /c [5] 

7 rV = ] n-x + 0(k , 34 

where y~^{x) G d^l is the preimage of n = (a; — k)/\x — k\ under the Gauss map, and 
^{y~^) is the Gauss curvature at G Q. Formula (l34l) has asymptotic character, and 
we want to figure out the range of values of k that give a good approximation of 7. We 
also analyze the approximation of the far field pattern f{k, x) by the measurement of 
the scattered field at the distance r. 

In order to conduct a numerical experiment, we restrict ourselves to the case where 
the direct problem can be easily solved. For that purpose we consider the problem that 
has an exact solution - scattering of plane wave e**^^ from a sphere of radius a with 
constant boundary impedance. Then the formula flM|l takes the form 



a + 2\f(k,x)\ .9 ^,,1, TT ^ 

where 6 is the polar angle on sphere. 

Similar to ([1]), we need to solve the boundary value problem for the Helmholtz 
equation 

— Am = k'^u, r > a, (36) 
Ur -ik'yu\^^^ = 0, (37) 

where 7 > is a constant surface impedance and 

u{x) = e'''' + ip{x) (38) 

with (f{x) satisfying the radiation condition (jl]). 
Solution of the problem is given by 

(39) 

where jn{z) and hn\z) are spherical Bessel functions of the first and third kind, respec- 
tively, and Pn{x) are Legendre polynomials [8]. From this formula we can determine the 
far field pattern f{k,x) (jl]) and calculate the surface impedance 7 using it asymptotics 
(!35|) for large values of k. 
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Figure 1: Reconstruction of the boundary impedance 7 = 2 of the illuminated part of 
the sphere of radius a = 1 from the formula ( l35l) . Left: the far field pattern f{k,x) 
with the wave number k = 200 is approximated by the amplitude of the scattered 
wave at different distances r from the center of the sphere. Right: the far field pattern 
f{k,x) is approximated as before with r = 200a, wave numbers k vary. Deviation 
from 7 = 2 increases as angle 6 approached 90° where incident rays are tangent to the 
sphere. 



6 Numerical examples 

Figure [T] shows the reconstructed boundary impedance 7 = 2 of a unit sphere based 
on the asymptotic formula (l35l) . In the left figure, the far field pattern with the wave 
number k = 200 was determined from the exact solution using (jl]) 

\f{k,x)\^\x\^{x) (40) 

for various distances |a;| from the sphere. The accuracy of approximation monotonically 
improves as the polar angle 6 increases and does not exceed about 0.5% for the distances 
beyond r = 100a. The right figure shows the dependence of the restored impedance 
on the wave number of the incident plane wave while at the distance r = 200a from 
the sphere. As the wave number k decreases, not only approximation of 7 deteriorates, 
but it also starts to exhibit oscillatory behavior. Approximation of 7 is also improving 
for larger values of 6 and remains below 0.5% as long as k > 50. 

The above approach leads to a good approximation of the impedance if the far 
field is measured at a distance by order of magnitude greater than the diameter of the 
sphere and for the wave length that is by order of magnitude lesser than the diameter of 
the sphere. Similar results are observed in restoring a compactly supported boundary 
impedance of a half space. Using the measurements of the far field pattern and both 
formulas fl33|) . f|T^ . and f|T2l) . we reconstructed the boundary impedance [2]|^a). In figure 
[21 we used explicit formula f|T^ with /c = 10 in (b). Then the Fourier coefficients of 
the far field pattern were perturbed by random numbers uniformly distributed in the 
interval [—1,1]. Figures [2](c)-(d) show reconstructed boundary impedance for A; = 15 
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when the amphtudes of the additive random noise were 1% and 5% of the greatest 
Fourier coefficient, respectively. 




Figure 2: Reconstruction of the boundary impedance 7(x) (a) using exphcit formula 
formula f|T^ for different values of the wave vector k. In (b) k = 10. k = 15 with 1% 
(c) and 5% (d) additive random noise. 

Reconstruction of the boundary impedance by asymptotic formula (!33|) gives slightly 
lesser accuracy as compared with exact formula (1141) . In figure E] we reconstructed the 
boundary impedance from figure using asymptotic formula fl33l) and k = 15. The 
Fourier coefficients of the far field pattern are corrupted by an additive uniformly dis- 
tributed random noise form the interval [—1, 1] with the amplitude 1% (a) and 5% (b) 
of the largest Fourier coefficient. The accuracy of reconstruction in this case is slightly 
less as compared with exact formula (|T4l) . 

Finally, in figure H] we reconstruct the impedance in the presence of a 1% additive 
random noise when the wavenumber is as small as = 5. Although in this case there 
is a significant error in the restored amplitude, it captures qualitatively the shape and 
the location of the inhomogeneity of the impedance. 
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Figure 3: Reconstruction of the same boundary impedance 7(x) from figure [2^ using 
asymptotic formula fl33l) when k = 15 and the amphtude of the additive random noise 
is 1% (a) and 5% (b), respectively. 

7 Conclusions 

We have considered the problem of determining a compactly supported boundary 
impedance from knowledge of the time harmonic incident wave and its far field pat- 
tern. The approach is based on a special selection of the directions in which the far 
field pattern is measured. Then the boundary impedance is expressed through a poten- 
tial using a simple exact formula, while the Fourier coefficients of the potential equal 
the measured far field patterns. Efficiency of the approach is illustrated by numerical 
examples. 
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Figure 4: Reconstructed boundary impedance 7(x) from figure [2^ using either explicit 
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